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Abstract 

We study exclusive Xc0,i,2 decays to four-hadron final states involving two charged and two 
neutral mesons: 7r + 7r~7r°7r°, K + K~7r°7r° , ppTT°7r°, K + K~7]ir° , and i^ ± 7r =F K°n° . The \c states are 
produced in radiative decays of 3.08 million if)(2S) resonance decays and observed in the CLEO 
detector. We also measure the largest substructure contributions to the modes 7r + 7r~7r°7r° and 
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Exclusive charmonium decays have been a subject of interest for decades as they are an 
excellent laboratory for studying quark-gluon dynamics at relatively low energies. However, 
current measurements in the P-wave Xc sector are sparse pj. Although these states are not 
directly produced in e + e~ collisions, they are copiously produced in the radiative decays 
ip{2S) — > 7Xc, each of which has a branching ratio of around 9% Recent data taken by 
the CLEO detector at the Cornell Electron Storage Ring to study e + e~ annihilations with 
a center of mass energy corresponding to the ip(2S) mass allow for a detailed study of Xc 
decays. 

Past research indicates that the Color Octet Mechanism (COM) plays a role in the decay 
of these P-wave charmonia states [1, 0, 0, @, 0] . In order to build a comprehensive under- 
standing about the P-wave dynamics, both theoretical predictions employing the COM and 
new precise experimental measurements for Xc many-body final states are required. Fur- 
thermore, decays of Xc, in particular Xco,2, may provide a window on glueball dynamics |8|. 



This analysis follows the general method of our earlier work on rj, rf [9J and three-body [l0(] 
decays of Xc, extending it to higher multiplicity states. The four-body exclusive Xc decay 
modes studied in this article contain two neutral and two charged hadrons in the final state 
and are being measured for the first time. We also take a first look at the gross features of 
the rich substructure in these decay modes. 

The data used in this analysis consist of 2.74 pb" 1 and 2.89 pb" 1 , a total of (3.08±0.09) 
xlO 6 ip(2S) decays, taken with the CLEO III [llj and CLEO-c [12| detector configurations, 
respectively. Both detector configurations provide 93% solid angle coverage, and have com- 
mon components of particle identification which are critical to this analysis: the main drift 
chamber, the Ring Imaging Cherenkov Detector (RICH), and the Csl calorimeter (CC). We 
distinguish two regions of the CC in polar angle: the barrel (| cos(0)| < 0.81) and the end- 
cap (| cos(#)| > 0.81). The CC detects photons with an energy resolution of 2.2% (5%) for 
photons with energy of 1 GeV (100 MeV), making it possible to resolve the three Xc states. 

The event reconstruction and final state selection criteria proceed along the lines of 
our recent analyses on this subject 0, [h|. We reconstruct the following Xc decay modes: 
K + K~r]7[°; K ± 7r Zf KgTT ; and h + h~7T 7r°, where h = n,K,p. In addition, we reconstruct the 
transition photon from ip(2S), thus detecting the entire event. 

Each charged particle in the event is required to pass standard criteria for track quality 
and geometric acceptance. We also require the number of such tracks to be either two or 
four depending on the final state. We demand all tracks come from the beam spot with 
a momentum-dependent cut on the impact parameter. This cut is less restrictive for low- 
momentum tracks, for which the resolution is poorer. We combine ionization loss in the 
drift chamber (dE/dx) and RICH information to discriminate between p, K, and 7r using 
X 2 criteria discussed elsewhere (l3j . Additional requirements suppress charged lepton QED 
backgrounds. We reject electron candidates as follows: for all tracks, we compute the ratio of 
CC energy to track momentum, Eqq/p, and the difference between the measured dE/dx and 
the expected dE/dx for the electron hypothesis, normalized to its standard deviation, <j e . 
We reject tracks with both 0.92 < E cc /p < 1.05 and \a e \ < 3. Particles that penetrate 
more than five nuclear interaction lengths of the muon detectors are rejected. The particle 
identification together with the lepton veto criteria are found to be more than 97% efficient 
for all the modes. 

We define photon showers as those having a lateral profile in the CC consistent with 
a photon and possessing at least 30 MeV of energy. We require photon candidates found 
in the endcap CC region that are used in tt° and rj reconstruction to possess more than 
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50 MeV of energy. We reconstruct 7r° — > 77 and 77 — > 77 candidates using a pair of photons 
that are kinematically fit to the nominal 7r° or 77 mass using the event's charged tracks to 
define the origin of the photon trajectories. A cut is placed on the mass fit of x 2 < 10 
(for one degree of freedom). Each tt° and 77 photon daughter is forbidden to be part of 
any other final state particle. We also reconstruct rj — > 7r + 7r~7r° decay by combining a pair 
of charged pion candidates with a 7r° and doing a mass-constrained fit to the nominal 77 
mass with the requirement x 2 < 10 f° r one degree of freedom. We reconstruct the Kg by 
constraining a pair of oppositely charged pions to come from a common vertex. We require 
that the reconstructed invariant mass of the two pions be within 10 MeV (~ 3.2 standard 
deviations) of the nominal Kg mass. 

We combine any unused photon with the four hadrons to reconstruct the complete event 
which we then constrain to the four-momentum of the ip(2S), using the nominal mass of the 
ip(2S) and taking into account the crossing angle between the e + and e~ beams (~ 4 mrad). 
We demand x 2 < 25 (for four degrees of freedom) for this constraint. This requirement 
strongly rejects background, and the fitting procedure greatly improves the mass resolution 
of the Xc- In < 10% of the events (the 77r + 7r~7r°7r° state is worst), multiple possible pairings 
of photons lead to multiple Xc candidates. We choose the one with the smallest x 2 ■ 

We study the efficiencies and resolutions of the final states by generating signal events 
using a GEANT-based [13] detector simulation. Over 0.7 million events were generated 
for the two detector configurations, three Xc mesons, and five final states. Events were 
generated in accordance with an electric dipole (El) transition production cross section of 
1 + A cos 2 9, where 9 is the radiated photon angle relative to the positron beam direction, 
and A = 1, —1/3, +1/13 for J = 0, 1, 2 particles, respectively. 

The efficiencies averaged over the CLEO III and CLEO-c datasets (weighted by the num- 
ber of ip(2S) events) are listed in Table [J for each final state. The efficiency includes the 
V 77 > V 7r + 7r~7T° and Kg — ► 7r + 7r~ branching ratios [lj]. The invariant mass distri- 
butions of the final state hadrons were fitted to three signal shape functions corresponding 
to each of the three Xc states and an additional constant background function. Each signal 
shape function consisted of a Breit-Wigner convolved with a double Gaussian resolution 
function. The widths of the Breit-Wigner functions are fixed to the intrinsic widths of the 
Xc states: T Xc0 = 10.4 MeV, T Xcl = 0.89 MeV, and T Xc2 = 2.06 MeV Q. The detector 
resolution is obtained from simulation by a fit to the difference between the generated and 
reconstructed invariant mass of the Xc products, mode-wise and separately for each Xc- The 
detector resolution, ranging from 4.5 to 7.1 MeV for the various final states, is less than the 
natural width of the Xco, but greater than the natural width of the xa and Xc2- The Xc 
masses are fixed to their nominal values [l| during the fitting. In all cases, the reconstructed 
masses, when allowed to float, are consistent with the expected values. 

Clean signals of Xco, Xd an d Xc2 are found in most of the modes studied as seen in 
Figure [TJ The signal yields and efficiencies, e, are listed in Table [I] for each final state, 
calculated assuming four-body phase space. The errors on the efficiencies due to limited 
simulation statistics are considered with other systematic uncertainties listed in Table IIHI 
The yields and efficiencies in Table |T] are used to calculate the final branching fractions for 
all modes, except 7r + 7r~7r°7r° and K + 7!-~K°tt°, which we proceed to study in more detail. 

We can further investigate how the decays proceed by searching for substructure in the 
four-hadron final states. Substructure can also affect the detection efficiency. Here we re- 
strict ourselves to looking for the gross features of substructure by plotting the invariant 
mass of di-hadron combinations in the Xc signal regions, after subtracting the yields from Xc 
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FIG. 1: The fitted distributions of the invariant mass of Xc candidates in data. The fitting procedure 
is described in the text. 



sidebands. The signal regions for the invariant mass distributions of Xco, Xci, and Xc2 can- 
didates, are 3.370-3.470 GeV/c 2 , 3.490-3.535 GeV/c 2 and 3.535-3.590 GeV/c 2 , respectively. 
One-sided invariant mass distribution sidebands of 3.300-3.360 GeV/c 2 , 3.470-3.490 GeV/c 2 
and 3.590-3.625 GeV/c 2 , respectively, are subtracted. We claim and show signals only for 
those cases where we observe a signal for resonant substructure with significance of > 4 a; 
these are all in the four body-modes 7r + 7r~7r 7r° and K ± n Zf K 7r°. After identifying the in- 
termediate states, we fit the distributions with Breit-Wigner functions for each of these 
intermediate states using a fixed mass and width as listed in and add an additional 
third order polynomial function to describe the background. We obtain the efficiencies for 
detecting the intermediate states, e', by fitting the resonant signals using the procedure de- 
scribed above, in samples of events that simulated the decay as proceeding entirely through 
the respective intermediate states. Figures [2] and [3] show two-body intermediate states in 
7r + 7r~7r°7r° and K n T KgTT° final states, respectively, for each of the Xc states. In the for- 
mer case there are two combinations entering the plot for each \ c candidate. We find clear 
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TABLE I: Yields (N) and efficiencies e for four-hadron final states. Efficiencies were estimated 
using a Monte Carlo sample generated according to four-body phase space as described in the text. 



Mode 


XcO 


Xcl 


Xc2 




N e(%) 


N e(%) 


N e(%) 


7r + 7T - 7T 7r 

K + K-ir°ir 

pp-K°ir 

K+K-r/TT 


1751.4±51.3 17.6 
213.5±16.8 12.7 
39.5±8.5 12.8 
56.4±9.2 6.19 
401.7±22.4 10.3 


604.7±28.5 17.5 
45.1±8.5 13.1 
11.5±4.2 13.9 
21.0±5.7 6.17 

141.3±13.3 11.2 


903.5±32.6 17.6 
76.9±9.7 12.6 
29.2±5.9 13.3 
22.9±6.1 6.07 

211.6±15.4 10.3 



signals for the p and K* mesons. The yields and efficiencies, e', of final states including 
these resonances are listed in Table [III and are used as inputs to the final branching fraction 
calculations. 

Several sources of systematic uncertainty in the branching fractions are listed in Ta- 
ble [TTTJ These include uncertainties in the number of ip{2S) particles, determined according 
to the method described in 0], tracking and particle identification efficiencies associated 
with charged particles, reconstruction efficiencies due to simulation statistics, and trigger 
simulation. By comparing the effect of varying the requirement on the \ 2 of the constrained 
four-momentum fit in data and in simulated events, we estimate a ±4.0% systematic error 
in modelling the \ 2 distribution. The uncertainty in the efficiency of the transition photon 
reconstruction for each of the final states is taken to be 2%. In addition, the systematic 
uncertainties for 7r° and rj efficiencies are 4% for each tc° or rj meson in the final state. Kg 
reconstruction introduces an additional uncertainty of 2% in one of the modes. The ro- 
bustness of the fitting procedure was checked by systematically re-fitting the Xc invariant 
mass plots using 1 o variations of masses, widths, and resolutions. The uncertainties esti- 
mated from this study vary between 0.2% and 0.7%. To account for the systematic error 
arising from deviations in the angular distributions of final states due to the presence of 
intermediate resonances, a sample of events that simulated the decay as proceeding through 
the respective intermediate states, was generated for the substructure modes: /o(980)7r°7r°, 
/ 2 (1270)7r vr , p ± 7T T 7r° and / (980)7r+7r- for the non-resonant mode 7r+7r-7r 7r°; K^K^tt 
and f (980)K + K- for the non-resonant mode K+K'tt ^ ; K*°K s n°, K^K^, K^K^n , 
p ± K^K° s , K*°K ± tc t and ^(1270)^° for the non-resonant mode K^K^w . We obtained 
the efficiencies of both the resonant and non-resonant modes by fitting the \ c signals obtained 
from MC events with and without substructure. Based on the differences in efficiencies be- 
tween final states with and without intermediate resonances, and assuming that there can 
be no additional unobserved resonances that can be more than 50% of the signal, we esti- 
mate the systematic error to be between 3.5% and 7.4% for the modes we studied. For the 
remaining modes, we conservatively assign a 7.5% systematic error based on an assumption 
that up to 75% of our events may contain substructure, and that the difference in efficiency 
of these resonant events to the non-resonant is no more than 10%. 

Systematic uncertainties for the intermediate state resonance modes include common 
sources to their non-resonant final states. Uncertainties were added in quadrature, and 
those due to the branching fractions of ip(2S) — > 7% c |2| are quoted separately. We use the 
values of B(ip{2S) — > 7% c ) in [2|, as those measurements use a subset of our data and similar 
analysis criteria, thereby enhancing the cancellation of systematic errors. For all final states 
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FIG. 2: Intermediate states for the decay Xcj — ► vr + 7r tt°tt°. Invariant mass combinations of (a) 
tt+tt for J = 0, (b) vr-yr for J = 0, (c) vr+vr for J = 1, (d) vi-vr for J = 1, (e) vr+vr for J = 2, 
and (f) vr _ vr° for J = 2 are shown after sideband subtraction. All plots show significant evidence 
for production of p^. 



except 7T + 7r 7T°7r and K ± tc t K s tt°, we convert the yields in Tables [J and [TT] to branching 
fractions using: 



Ni 



cJ 



(1) 



where Ni is the yield; N^s) is the number of ip{2S) = 3.08 xlO 6 ; e is the same e listed 
in Table [TJ for four-hadron modes, or e' listed in Table [UJ for substructure modes, and i 
represents a particular decay mode. 

To calculate the branching fractions for the inclusive four-hadron final states for 
7r + 7r - 7r°7r° and K ± Tr Zf K°7i° we use a modified procedure. Since the p ± ii T 7i clearly domi- 
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FIG. 3: Figure shows intermediate states for the decay Xcj —* K^tt^ Kgir . Invariant mass combi- 
nations of (a) K^ir? for J = 0, (b) K^ir? for J = 1, (c) K^ir* for J = 2, (d) v^vr for J = 0, (e) 
tt ± tt for J = 1, (f) t^tt for J = 2, (g) K±tt° for J = 0, (h) ^tt for J = 2, (i) Kgvr for J = 2 
(j) Tr^/fg for J = 0, and (k) ir^Kg for J = 2 are shown after sideband subtraction. Strong signals 
for or K* production are visible. 

nates the four-hadron final state yields for the 7r + 7r~7r°7r° mode (Tables HI and HT|) . we use 
the efficiency s of the p ± n T 7i sub-mode listed in Table HT1 and Equation [TJ to determine the 
7r + 7r~7r 7r° branching fraction. The efficiencies e were obtained by fitting the Xcj signals in 
substructure simulations using the same fitting procedure as that used for the four-hadron 
signal simulations 1 . To calculate the K^tt^K ^ branching fraction, we modify the proce- 
dure by taking into account that this channel has many intermediate resonances (Table HT1). 
We replace the ratio iVj/ej in Equation [TJ with an efficiency corrected yield (Y K ± 7r ^ K o 7r o) 
by adding the individual efficiency corrected contributions of all resonant and non-resonant 



The efficiency e' is lower than e since the sideband subtraction procedure used for obtaining the efficiency 
e' results in some loss of efficiency. 
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TABLE II: Yields and efficiencies (in %) for substructure modes, e represents the efficiency obtained 
by fitting the \ c signals using the same procedure as that used for the non-resonant four-hadron 
modes, in substructure simulated samples. Efficiency e' was obtained by fitting the intermediate 
peaks after applying the sideband subtraction procedure described in the text. N' describes the 
yield corresponding to this efficiency e' . 



Mode 


XcO 


Xcl 


Xc2 




N' e e' 


N' e e' 


N' e ef 


p~7r + 7r° 


661.4±56.8 17.4 15.7 
697.1±56.6 17.4 15.7 


355.7±42.1 16.9 16.3 
356.6±41.0 16.9 16.3 


519.3±38.7 16.8 16.1 
512.6±40.0 16.8 16.1 


K* K° s ir°, K*° K±ttT 
p ± K^K° s 


52.5±12.1 11.0 9.94 

64.1±12.8 10.3 9.30 
42.5±10.0 10.4 9.44 
179.7±22.7 11.0 9.92 


37.9±8.5 11.1 10.9 
79.5±16.9 10.8 10.6 


63.0±10.5 11.5 11.2 
38.7±9.0 9.50 9.14 
51.1±9.8 9.66 9.40 
39.3±8.7 9.44 9.10 
62.9±15.9 10.8 10.5 



TABLE III: Systematic uncertainties (fractional errors in %). The overall systematic uncertainty 
is obtained by adding the individual contributions in quadrature, with the exception of the photon 
simulation uncertainty which adds linearly with the 7r° and rj reconstruction uncertainties. 



Source 


7T+7T 7r°7r° 


K + K--K ir° 


pp7r°ir° 


K + K-r]7T° 




N i>{2S) 


3.0 


3.0 


3.0 


3.0 


3.0 


Tracking efficiency 


1.4 


1.4 


1.4 


1.4 


2.8 


Particle identification 


0.6 


2.6 


2.6 


2.6 


2.2 


Simulation statistics 


1.9 


2.3 


2.3 


2.8 


4.5 


Trigger efficiency 


1.0 


1.0 


1.0 


1.0 


1.0 


Kinematic constraint cut 


4.0 


4.0 


4.0 


4.0 


4.0 


Transition 7 simulation 


2.0 


2.0 


2.0 


2.0 


2.0 


7r°,7/ — > 77 reconstruction 


8.0 


8.0 


8.0 


8.0 


4.0 


Kg Vertexing 


0.0 


0.0 


0.0 


0.0 


2.0 


Fitting procedure 


0.5 


0.6 


0.7 


0.7 


0.2 


Model dependence 


3.5 


7.4 


7.5 


7.5 


6.0 


Overall 


12.0 


14.0 


14.0 


14.1 


11.6 



channels computed as: 

ik^kw = \- 2^ — r ' 1 i 2 ) 

^_fsr±7rTK 7r° k= i \ ZK+tt^ROtt / 

where k runs over the substructure modes we consider, K*°Kgii , K*° K^tc 11 , K^KT-k , 
K^tt^K®, and p ± K^K° s . 

The branching fractions obtained and their uncertainties are summarized in Table HVl 
Where we do not find evidence of a signal, we present a 90% CL. upper limit by determining 
the value that includes 90% of the probability density function (p.d.f ) obtained by convolving 
the p.d.f for the branching fraction with a Gaussian systematic error. 

The results of the branching fractions of the four-hadron modes studied in this analysis 
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include both resonant and non-resonant contributions. Our three-hadron intermediate res- 
onance branching fractions are inclusive and may contain additional resonant substructure 
that we do not explicitly measure; therefore, the branching fraction measurements cannot 
be trivially related to the amplitude for the specific three-body non-resonant decay. We 
assume variations of the efficiency due to this additional substructure produce changes in 
the branching fractions within the model dependence systematic error. 

We compare the branching fractions Xcj - > p ± tt t tt° for J = 0,1,2 (Table IIVD with 
measurements of Xcj — * p°tt + tt~ [if, and observe that our results are consistent with a ratio 
equal to unity as expected from isospin conservation. 

Furthermore, we also find our measurement of B(x C 2 K*°K ± n T ) (Table |V|) to be 
consistent within experimental errors with B(x C 2 K*°K ± ir T ) of [l|. Our results are con- 
sistent with the expected isospin-related prediction that \ c — > K*Kit where \ c — > K*°K°7r° 
and Xc K^K^n have equal partial widths. Moreover, our measurements for B{xc2 
K^KV), B(xco -> K^n^K ), and B{ Xc 2 -^K*VK°) (Table EJ) are in good agreement 
with the isospin expectation of 

B( X c -> K*°K\°) : B( X c -> K^K^T) = 1:2 (3) 

and 

B( X c -> K*°K°7T ) : B{ X c -> K^K ) = 1:2. (4) 



TABLE IV: Branching fractions B with statistical and systematic uncertainties are shown. The 
symbol "x" indicates product of branching fractions. The third error in each case is due to the 
ip(2S) — > 7Xc branching fractions. Upper limits shown are at the 90% C.L and include all the 
systematic errors. The first line of each set of measurements includes the contributions from the 
substructure listed below. 



Mode 


XcO Xcl 


Xc2 




B (%) B (%) 


B (%) 


7r + 7r~7r 7r 


3.54±0.10±0.43±0.18 1.28±0.06±0.15±0.08 


1.87±0.07±0.22±0.13 


p + 7T~7T 
p~7T + 7T 


1.48±0.13±0.18±0.08 0.78±0.09±0.09±0.05 
1.56±0.13±0.19±0.08 0.78±0.09±0.09±0.05 


1.12±0.08±0.13±0.08 
1.11±0.09±0.13±0.08 


K + K-ir°ir 


0.59±0.05±0.08±0.03 0.12±0.02±0.02±0.01 


0.21±0.03±0.03±0.01 


pp7t°7r° 


0.11±0.02±0.02±0.01 < 0.05 


0.08±0.02±0.01±0.01 


K+K-r/Tr 


0.32±0.05±0.05±0.02 0.12±0.03±0.02±0.01 0.13±0.04±0.02±0.01 


K ± TT T K°TT° 


2.64±0.15±0.31±0.14 0.92±0.09±0.11±0.06 


1.41±0.11±0.16±0.10 



K*°K°tt°x K*° - 
F^Vx K*° - 

p^K^K 



K*** 0.37±0.09±0.04±0.02 0.25±0.06±0.03±0.02 0.39±0.07±0.05±0.03 
+ K°*k° - - 0.30±0.07±0.04±0.02 

-» ir±K° 0.49±0.10±0.06±0.03 - 0.38±0.07±0.04±0.03 

-» l^vr 0.32±0.07±0.04±0.02 - 0.30±0.07±0.04±0.02 

1.28±0.16±0.15±0.07 0.54±0.11±0.06±0.03 0.42±0.11±0.05±0.03 



7T + 7T 7T 7T , 



In summary, the branching fractions for Xcj 
K + K-r]7r°, and X cj -> KVK !? for J = 0, 1, 2 and, Xc j 
sured for the first time. For the mode Xci 



XcJ 



K+K-n 7i°, Xcj 

0, 2 are mea- 



ppii TT for J 



pp7i n° upper limits at 90% C.L. are presented. 
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TABLE V: Branching fractions and total error measurements for the isospin-related K*Kir inter- 
mediate modes. 



Mode 


XcO 


Xcl 


Xc2 




B (%) 


B (%) 


B (%) 




0.56±0.15 0.38±0.11 0.59±0.14 








0.90±0.25 




0.74±0.18 




0.57±0.13 


K* ±r K T K® 


0.96±0.25 




0.90±0.25 



We also measure for the first time the partial branching fractions of Xcj vr + 7r°, Xcj —> 
K* K°n°, and X cj -> p ± K^K° for J = 0,1,2; X c2 -> K*° K^; Xcj -> i^* ± -f^ =F 7r°and 
Xcj — ► K* ± n T K° for J = 0, 2. These four-hadron final states account for up to 8% of the 
hadronic width of the Xc states. 

These measurements improve our existing knowledge of the exclusive multi-body decay 
modes of the Xc states and provide insight into their decay mechanisms. The four-hadron fi- 
nal states 7r + 7r~7r°7r° and K ± 7r^K°7r° contain a rich substructure of intermediate resonances. 
These interesting results form a basis from which higher statistics studies of hadronic Xc 
states and their substructure may follow. 
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